Qubit Coherent Control with Squeezed Light Fields 
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We study the use of squeezed light for qubit coherent control and compare it with the coherent 
state control field case. We calculate the entanglement between a short pulse of resonant squeezed 
light and a two-level atom in free space and the resulting operation error. We find that the squeezing 
phase, the phase of the light field and the atomic superposition phase, all determine whether atom- 
pulse mode entanglement and the gate error are enhanced or suppressed. However, when averaged 
over all possible qubit initial states, the gate error would not decrease by a practicably useful amount 
and would in fact increase in most cases. We discuss the possibility of measuring the increased gate 
error as a signature of the enhancement of entanglement by squeezing. 

PACS numbers: 42.50.Ct, 03.67.Bg, 03.65.Yz, 32.80.Qk 



I. INTRODUCTION 

The coherent control of a two-level atom in free space 
with resonant laser fields is important both from a fun- 
damental as well as a practical standpoint. In particu- 
lar, for quantum information processing purposes, where 
the two- level atom encodes a single quantum bit (qubit), 
coherent control accuracy plays a crucial role. Fault- 
tolerant quantum error correction schemes, necessary for 
large scale quantum computing, demand that the error in 
a single quantumgate is below a certain, typically very 
small, threshold [lj. 

In recent years the fundamental limitation to qubit 
control accuracy imposed by the quantum nature of the 
control light fields has been the focus of several studies 
@, 1,11 ESQ. For a coherent state of light (CL), where 
both phase and amplitude fluctuations are at the stan- 
dard quantum limit (SQL), it was found that the gate 
error is determined by the qubit spontaneous decay rate. 
At a time shorter than the qubit decay time, part of that 
error can be interpreted as due to entanglement between 
the control field mode and the qubit (forward scattering) . 

This entanglement is a measure of the information on 
the atomic state carried away by the pulse field and could, 
in principle, be exploited as a resource in quantum com- 
munication schemes. For a CL pulse of duration r, the 
entanglement, measured by the tangle T, was found to 
be T ~ r ~ where n is the average number of 

photons in the pulse. Thus, in agreement with the cor- 
respondence principle, quantum entanglement becomes 
very small as the pulse approaches the classical limit, i.e. 
n > 1. 

It is therefore natural to ask whether a control with non 
coherent light held states would reduce or increase both 
the operation error and atom- light mode entanglement. 
In this paper we investigate the fundamental quantum 
limitation to the fidelity of qubit coherent control with 
squeezed light (SL) fields. As opposed to CL, SL has 
phase-dependent quantum fluctuations, where one field 
quadrature has reduced and the other increased fluctu- 
ations compared with the SQL. Reduced phase fluctua- 



tions of squeezed light have been shown to, in principle, 
improve the fundamental limitation to phase metrology 
accuracy Q . The use of squeezed light states was also 
shown to be beneficial for continuous variable quantum 
key distribution [l(| El, EH- Here, we explore the pos- 
sible advantage of the use of squeezed light for qubit co- 
herent control. Both the error of a single-qubit quantum 
gate and the entanglement between the atomic qubit and 
the SL mode are calculated. 

Atom-SL interactions were studied extensively in the 
past [l3|, In most cases the atomic dynamics was 

found to be modified in a squeezing-phase-dependent 
way. Gardiner [l5| found that a two-level atom damped 
by a squeezed vacuum reservoir exhibits squeezing-phase- 
dependent suppression or enhancement of atomic coher- 
ence decay. Carmichael [n| added a coherent field that 
drives the atom and showed that the dynamics of the 
atomic coherence depends on the squeezed vacuum and 
the coherent driving field relative phase. Both assumed 
a broadband squeezed vacuum reservoir whereas the fi- 
nite bandwidth case was studied in [17]. Milburn, on 
the other hand, studied the interaction of a single-mode 
SL with an atom using the Janes-Cummings model and 
found a squeezing-phase-dependent increase or decrease 
of the collapse time [18| ■ 

Here we assume a qubit that is encoded in an atomic, 
two-level, superposition and a control field that is realized 
with a pulse of SL in a single electromagnetic (EM) field 
mode and of duration shorter than the typical atomic de- 
cay time, in free space. We find that the squeezing phase, 
the phase of the light field and the atomic superposition 
phase, all determine whether atom-pulse mode entangle- 
ment and the gate error are enhanced or suppressed. It 
is found that, although reduced for certain qubit initial 
states, when averaged over all possible states the minimal 
gate error is comparable to that with CL fields. 

The paper is outlined as follows. Section ILTl describes 
a general model of the interaction between a light pulse 
and an atom in free space within the paraxial approxi- 
mation. In Sec. IIIII the pulse mode is assumed to be in 
a squeezed state and its effect on the atom is examined. 



2 



Atom-pulse entanglement as well as the operation error 
are calculated. Some aspects of possible experimental re- 
alization are discussed in Sec. IIVI where our results are 
also briefly summarized. 



II. GENERAL MODEL AND APPROACH 

We describe the interaction of a two-level atom with a 
quantized pulse of light propagating in free space. Classi- 
cally the propagation of such a pulse is well approximated 
by a paraxial wave equation. A paraxial formalism for 
quantized pulses was described in [I El HI- We briefly 
describe theparaxial model following a similar approach 
to that in J4[. The derivation of this approximation is 
outlined in the appendix, where we start from the ba- 
sic atom-field Hamiltonian in free space and get to the 
model described below. 



A. The pulse mode 

Consider a laser pulse traveling in the z direction and 
interacting with a well-localized atom or ion at position 
f a ■ This pulse is a wave packet of Fourier modes around 
some wave vector ko = koz, where z is a unit vector 
pointing to the z direction, i.e. it is composed of a car- 
rier wave vector ko and a slowly varying envelope. The 
pulse is assumed to propagate with a small diffraction 
angle, so that the paraxial approximation is valid. A 
natural description for such a pulse is as a superposition 
of Gaussian beams with a well-defined transverse area 
A at the position of the atom, each with wave number 
kg + k, where k <C ko, analogous to a narrowband wave 
packet of Fourier modes in Id. In this approximation the 
pulse is effectively a one-dimensional beam with trans- 
verse area A similar to a Gaussian beam. The envelope 
of the pulse mode, ipf (z), is a function of width cr where 
t is the pulse duration and c the speed of light. The field 
operator assigned to this pulse mode is o/ . 

B. Interaction picture and the atom-pulse bipartite 
system 

We start with the full atom-field Hamiltonian in free 
space, 

H = H F +H A + Haf 

Ha + Hp = — hw a a z + fio^a^a^ 

k /( 

H AF =^ih(g %iJ a % ^-h.c)) {a + + a_), (1) 

kfi 

where k[i are the indices of a Fourier mode and its po- 
larization respectively, and are the atom-mode cou- 



plings. Here the a operators are the spin 1/2 Pauli ma- 
trices and hui a is the energy separation between the two 
atomic levels. Moving to the interaction picture with 
respect to the interaction Haf in the rotating wave ap- 
proximation, assuming no detuning uo — koc — oj ai and 
using the paraxial approximation, we get the Hamilto- 
nian in the interaction picture (see appendix) 

Hi = Hs + Hsr 

Hs = ihgf (t)(a fo a + -a\ o a-) 

H S r = ih^2(g r ara + - g*ala-) , (2) 

r 

where the indices {r} denote all the free space EM modes 
apart from the fo pulse mode, with couplings g r which 
are in general time dependent. The fo mode coupling is 

.9/0 (*) = ^^A^Vfa (*« " (3) 

where z a is the atom's position. In Eq. © we assume 
both the envelope of the pulse <fif {z) and the dipole 
matrix element d to be real and set z a = 0. A similar 
model Hamiltonian was also used in 0, @| where the 
entanglement between a CL pulse and an atom was 
calculated. Note that the time dependence of the 
atom-pulse coupling <?/ (i) is identical to that of the 
pulse, the pulse therefore switches the interaction on 
and off. 

The relevant bipartite system is composed of the atom 
and the fo light mode, coupled by Hs- The rest of the 
EM modes except for fo, {r}, make up a reservoir that 
interacts with the bipartite system via the coupling with 
the atom. 



III. SQUEEZED LIGHT PULSE OPERATING 
ON AN ATOM 

We consider a pulse of squeezed light operating on an 
atom. The pulse's average field changes the atomic su- 
perposition coherently. In a Bloch sphere representation 
this corresponds to a rotation of the Bloch vector. The 
quantum fluctuations of the field lead to atom-pulse en- 
tanglement and add a random component to the average 
operation. We wish to analyze and calculate the error 
caused by quantum fluctuations as well as the degree of 
atom-pulse mode entanglement. The model outlined in 
Sec. II is sufficiently general to describe a pulse in any 
quantum state and its interaction with the atom. In 0, 0] 
calculations for CL state were performed. In this Section 
we will discuss quantum pulses of SL interacting with a 
two-level atom. The general formalism is presented in 
Sec. Ill A. In Sec. Ill B we consider a rotation of the 
Bloch vector by 180 degrees about an axis lying in the 
Bloch sphere equatorial plane (a n pulse) as an exam- 
ple and neglect the reservoir (i.e. the rest of the EM 
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modes). This allows to obtain analytic expressions for 
the atom-pulse entanglement which, in this case, is the 
sole generator of gate error. For longer pulses, reservoir 
effects become increasingly dominant for both entangle- 
ment and error generation. In Sec. Ill C we obtain the 
results in the presence of the reservoir, i.e. including 
spontaneous emission into all modes. 

A. Unitary transformations and master equation 

In order to determine the final state of the bipartite 
photon pulse-atom system in the presence of the reservoir 
of free EM modes, an appropriate master equation should 
be used. 



Here a = d/ , s = sinh(r) and c = cosh(r). In this 
representation the initial state of all the EM modes is 
the vacuum state. 

The Hamiltonian @ includes two types of driving 
terms. The second term in describes an atom driven 
by the average field and therefore contains no EM field 
operators. The first term in @ describes interaction with 
the EM squeezed vacuum. Unlike the coherent state case, 
it includes terms like d^cr+ similar to those that are typi- 
cally neglected in the rotating wave approximation. Here 
however, such terms appear as a result of squeezing, and 
indeed vanish for r = 0. These terms reflect the fact 
that unlike the ordinary vacuum the squeezed vacuum 
contains photons that can be absorbed by the atom. 



1. Initial state 

Here we assume the light mode /o to be a single pulse 
mode in a squeezed state and choose the following initial 
state for the EM field 



Here 



where 



W,e) fQ <g> |0) {r} . 
\a,e) fo =D(a)S(e)\0) fo , 



S(e) = S fQ (e) = e 



f - ea V 



(4) 
(5) 
(6) 



is the squeezing operator with the squeezing parameter 
e = re l2 ^ and 



D(a) = D fo (a) = e 



(7) 



is the displacement operator with the field average a = 
\a\e w . 

The phases 4> and 6 together with the phases associated 
with the initial state of the atomic qubit will play an 
important role in the following. 



2. Unitary transformations 

It is convenient to transform the problem into a repre- 
sentation where the initial state of the pulse mode is the 
vacuum state. To achieve this, we first use the Mollow 
transformation [2l|, [22j and move the average field a from 
the initial state into the Hamiltonian. This affects only 
the term Hs in the Hamiltonian ([2]) which becomes 

Hs = ihgf a (t) [a fa a + - h.c.]+ihg fo (t) [aa+ - h.c] . (8) 

We then make a time-independent Bogoliubov transfor- 
mation, S'(s)HS(£). This again affects only Hs which 
becomes, 



Hs = ihgf (t) [(7+ (ca — e 2l *sd^) — h.c] 
+ i h 9fo(t) l a °+ -a* a-]. 



3. Master equation 

We now examine a particular case of a n pulse of du- 
ration t. In a Bloch sphere representation and within 
our notations, the average field of such a pulse rotates 
the Bloch vector by 180 degrees about an axis, lying in 
the sphere equatorial plane, in an angle 8 from the y 
axis. For a constant interaction coefficient g this requires 
\a\gr = ir/2. The time-dependent interaction coefficient 
gf (t) from Eq. includes the pulse profile, which is 
assumed to be approximately constant during the time 
interval r and vanishing elsewhere. We set t = to be 
the time at which the pulse reaches the atom so at t = r 
it has already passed it, 



g fo (0<t<T)=g = 
g fo (t<0,t>T)=Q. 



(10) 



Here k = d 2 uJo/(2eohAc) is the atomic decay rate into 
all the paraxial modes of the same transverse profile as 
the fo pulse mode The dynamics during the pulse is 
then approximately governed by, 



H 

Hsr 



= H s + H, 



SI! 



= ihJ-cW + ih^- (e l9 a + - 
V t 2t 



(11) 



where W is the dimensionless operator that describes the 
quantum field part of the atom-pulse interaction, 



W 



cr_|_ I a — e 



a_ idt - e 2 ^-a 



(121 



The corresponding master equation for the bipartite sys- 
tem density matrix, p(t) , in the presence of the reservoir 
is 

^ = -^[Hs,p] - ^{v+a_,p} + r&-pa + , (13) 



(9) where T is the atomic decay rate in free space. 
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This master equation assumes a Markovian reservoir, 
i.e. a reservoir with a very short correlation time. Such 
assumption is not obvious for the free space EM reservoir 
that lacks a single pulse mode, like the one we consider. 
This rather artificial reservoir has a correlation time of 
the order of the pulse duration r. In section IV B of @, 
a similar Markovian equation was written for a CL pulse 
and was then solved perturbatively for kt <C 1. In fact, 
in [23| it was found that the assumption (gr) 2 =kt<1 
is actually a necessary condition for the validity of the 
Markovian assumption and the Master equation (TlU|) , to- 
gether with the condition k -C T. An intuitive argument 
for the requirement (gr) 2 — kt <^_1 can be based on the 
time scales related to the pulse mode that is excluded 
from the reservoir. While the correlation time of the 
reservoir is r, this correlation's effect on the dynamics 
should be of the order of n. Therefore Markovian dy- 
namics is obtained by requiring rer <C 1. Markovian time 
evolution guarantees that the probability of successive 
reabsorptions and reemissions of photons from and into 
the /o mode, typical of the Janes-Cummings dynamics, 
are negligible as it should be in free space. 

The initial bipartite state, p(Q) — \a, 0)(<r, 0|, where \a) 
denotes any atomic qubit state and /o is in the vacuum 0, 
is pure and separable. Solving Eq. (fi"3|) for t = t yields 
the bipartite density operator p(r) from which properties 
like atom-pulse entanglement can be calculated. 



B. Entanglement and error probability calculations 
neglecting the reservoir 

We start by considering a pulse duration which is much 
shorter than the atomic decay time An analytic so- 
lution to Eq. (TT5|) can be found in this case by neglecting 
the reservoir. In this short pulse limit, the probability 
to absorb a photon from the pulse and then to emit one 
to the reservoir is negligible. Then, one can indeed ne- 
glect any modification to the atom-pulse entanglement 
due to the presence of the reservoir. However, gate error 
calculated in this approximation is induced only by the 
quantum fluctuations of the pulse mode. We note that 
at least for the CL case, this error is a relatively small 
part of the total error which is mostly due to spontaneous 
photon emission to the reservoir EM modes. 



1. Bipartite wave function calculation 

Neglecting the reservoir, Eq. (fT5|) becomes a 
Schrodinger equation, the solution of which at t = r is 
|^( T )) = U(t)\o;Q), where 



U(t) = e 



(14) 



is the propagator for t = r. Here — is chosen and 
A = cosh(r) y/nr. We now take the assumption kt < 1 
one step further and require A < 1. This puts a limit 



on the magnitude of the squeezing parameter r. As will 
be shown later, in the paraxial approximation where k 
is small and with presently available squeezing technol- 
ogy, this is a very good approximation. Looking at the 
propagator in Eq. (|14p one notes that initially the ma- 
trix elements of W are O(l) so U{t) can be expanded in 
powers of A 



U(t) = fjW + Al>« + X 2 U^ + 



(15) 



where 



= l>o(£ = l) 



(16) 



(i 



rf6 / d£iU (l - a)WoKi - &)Wo(6) 



and where we defined Uo(0 = ex P [( 7r /2)(o'+ — ct_)£]. 

Using this approach we calculate the final states for 
various initial states of the atom. We find that up to 
0(A 2 ) it is enough to work in the truncated bipartite ba- 
sis that includes only 4 states: {\g, 0), \g, 1), |e, 0), |e, 1)}. 
States of the form |er, 2) do appear in the expansion. How- 
ever, in our calculations of the entanglement and gate 
error they do not contribute up to and including 0(A 2 ) 
and are therefore omitted. We thus end up with a trun- 
cated basis that, following the perturbative assumption, 
includes only a single photonic excitation in the pulse 
mode, as argued in Q. The resulting final state is 
not normalized, we therefore expand the normalization 
factor n(A) = 1/ \J to second order in A and use it 
to normalize the state. 



2. Entanglement calculation 

We measure the entanglement of the bipartite state 
using the tangle monotone which is related to the square 
of the concurrence 0, Hi| and compare our results for 
SL pulses with those obtained in 0] for CL pulses. 

For a pure bipartite state \ip) of two subsystems A and 
B and where one of the subsystems is a two-level system, 
the tangle takes the form [26j 



T(M(V>|) = 2(l-T^ [{Tt B U)(i>\)) 



(17) 



We calculate the atom-pulse tangle for various atomic 
initial states up to 0(A 2 ), by inserting the normalized 
final state, calculated using Eq. ([To]), into Eq. (TIT)) . 
For atomic initial states \g) and |e) we get 



1 



- 2 



cos(20) 



(18) 



For a CL pulse, i.e. s — and c = 1, the tangle is given 
by T m kt, which is the same as calculated in Q. For 
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non-zero squeezing this result is rather strongly modified 
and in particular the tangle depends on the squeezing 
phase 4>- As an example when e r 3> 1 and therefore s/c 
is close to unity, we obtain T s» 2c 2 kt[1 — cos(20)]. This 
expression has a typical interference pattern which will 
be discussed in section III B4. 

Let us now consider two different squeezing phases, 
= and 4> = it/2, which for 9 = correspond to 
amplitude and phase squeezing respectively. We obtain 
from (HU) 

T k Kre T2r , (19) 

where the minus (plus) sign is for <f> = (<f> = tt/2). When 
compared with the CL case, in the <fi = tt/2 case we get 
enhancement whereas in the <f> — case we get suppres- 
sion of entanglement. This result is plotted in Fig. [TJ 
where we set r = 0.1T -1 and k — T/1000. This estimate 
for k results from assuming a Gaussian beam focused on 
the atom with a lens of numerical aperture of about 0.1, 
consistent with the paraxial approximation 0, @| • With 
these numbers, the assumption A <C 1 is valid for r < 3. 
Current state-of-the-art light squeezing experiments do 
not exceed r = 1.16 (see Sec. IV), therefore A < 1 is a 
very good assumption. 

We repeat the tangle calculation for initial superpo- 
sition states of the form (e l8A \g) + |e)). In a Bloch 
sphere representation this state is represented by a vec- 
tor lying in the equatorial xy plane with an angle 9 a from 
the positive x direction. The tangle we get is 

p -i29 A 

T « [tt (1 + e a ° A ) e^ r + Ae l ° A e ±r ] kt, (20) 

where the upper (lower) sign refers to <p — (4> — tt/2). 
This expression is plotted in Fig. 1 for the two initial 
states 9a = and 9a = tt/2. For each initial state, 
a squeezing-phase-dependent enhancement and suppres- 
sion compared with the CL case are observed. 

As could be expected from atom-SL interaction, we 
observe dependence on the relative phase between the 
squeezing phase, 4», and the phase of the average of the 
pulse, 9 (which is chosen here to be 0). Furthermore, the 
tangle indeed depends on the atomic initial state. For 
fixed pulse phase, the change in entanglement for initial 
states that lie in the Bloch sphere equatorial xy plane 
depends on the difference between the squeezing phase 
and 9a- 

3. Error probability calculation 

The error probability of a quantum operation is just 
one minus the probability of the atomic qubit to be in 
the target state \rj) of the operation after the operation 
is performed. A classical light pulse can be chosen to 
transform an atom from the initial to a prescribed final 
state 1 77). When operating with the quantized pulse, the 




squeezing (r) 

FIG. 1: The tangle, Eq. (|17|) . calculated for the atom and the 
squeezed light pulse mode bipartite system vs the squeezing 
strength. The pulse duration is r = O.IF -1 and the coupling 
to the rest of the EM modes is neglected. The decay rate into 
paraxial modes is taken to be k = T/1000. Results are shown 
for four atomic initial states and for the squeezing phases 
cj> = 0,7r/2, corresponding to amplitude and phase squeezing 
respectively. Compared with the coherent light case (r = 
0), one obtains a squeezing-phase-dependent suppression and 
enhancement of the entanglement as squeezing gets larger. 



operation is subject to quantum fluctuations. The fidelity 
F is the probability of the atom to be in \rf) after the 
quantized field has operated on it, F = (t]\pa\t)), where 
Pa is the density matrix that describes the, generally 
mixed, final state of the two-level atomic system. The 
error probability is then, 

P error = 1 -F= 1 - (v\pa\v)- (21) 

Since in the absence of the reservoir the sole error is that 
due to atom-pulse entanglement, we were able to find a 
simple analytic relation between the tangle and the error 
probability calculated up to 0(A 2 ), T » 4P error . This 
linear relation implies that the results for P error follow 
similar trends to those for the tangle and are squeez- 
ing phase and initial atomic state dependent (Fig. 2). 
For quantum information processing purposes the rele- 
vant error is that averaged over all possible atomic ini- 
tial states, < P e rror >■ The averaged error is calculated 
analytically using the fidelity of 6 different initial atomic 
states [27| 

< P error >k (0.0675e ±2r + 0.1665e =F2r ) kt, (22) 

where as before the upper (lower) sign is for = 
((f) = tt/2). This result is plotted as a function of 
squeezing in Fig. 3. As can be seen the averaged er- 
ror, < Perror >i depends on the squeezing phase </>. For 
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1 1.5 2 

squeezing (r) 




0.5 1 1.5 2 2.5 3 

squeezing (r) 



FIG. 2: The error probability of an operation with a squeezed 
light 7r pulse on an atom vs squeezing strength. The pulse 
duration is r = 0.1T -1 and the rest of the EM modes are 
disregarded. As in Fig. 1, k = F/1000 is taken. The in- 
set shows the atomic initial states and the squeezing phase. 
Here <j> — (7r/2) corresponds to reduced amplitude (phase) 
fluctuations. The dependence on the squeezing phase follows 
the results for the entanglement. Lower errors than for the 
coherent light operation (r = 0) can be achieved. 



most values of the squeezing parameter r the average er- 
ror is larger than in the CL case (r = 0). This implies 
that on average, more entanglement is generated between 
the atom and the pulse in SL than in CL case. Still one 
observes that for the region < r < 0.45, and (f> = 
(amplitude squeezing), the average error drops relative 
to the error in the CL case. Though practically small, in 
this range of the squeezing parameter the SL operation 
has lower quantum limit to the operation error than with 
the CL pulse due to reduced intensity noise. 



4- Discussion of the results 

4a. Quantum interference 
Consider the interaction W between the atom and the 
fluctuating part of the light pulse shown in Eq. (fT2|) . 
In addition to the Jaynes-Cummings like atom-photon 
interaction term W\ = i{acr + — /i.e.), present also in 
the CL case, it contains a squeezing-dependent term 
W% = i [{s/c)e^ 2ip d- + a i — h.c.]. Atomic and light mode 

excitations can be generated via both W\ and W% . These 
two paths interfere quantum mechanically resulting a de- 
pendence of atom-light entanglement and error probabil- 
ity on both squeezing and atomic phases. 

Starting from \g, 0) we examine the bipartite state evo- 
lution under the first two operators in Eq. (|16p . The 



FIG. 3: Same as Fig. 2, but for the average error probability. 
The average is over all possible atomic initial states. For most 
cases, the average error with squeezed light increases relative 
to coherent light (r = 0). For the squeezing phase <f> — 
(amplitude squeezing) and squeezing strength < r < 0.45 a 
very small decrease in the average error is observed, as shown 
in the figure inset. 



classical field term, fj(°\ is evolving the state to |e,0) 
via superpositions of \g, 0) and |e, 0). Simultaneously the 
state evolves under the term W in U^. The W\ part of 
it couples only to |e, 0) evolving it to \g, 1). The W2 part 
on the other hand couples only to \g, 0) and transform 
it to |e, 1). Subsequent evolution under U forms super- 
positions of \e, 1) and \g, 1), the amplitudes of which are 
coherent sums (and therefore show interference) of evolu- 
tion amplitudes under both W\ and W^- Quantitatively, 



15,0} 



|e,0) + A / dtu (i-0(Wi+w 2 )u (0\g,o) = 



A 



,0)-Xj o d£l>o(l-£)sin(^) \g,l) 
J ^o(l-Oe- i20 ~cos(|£) |e,l) = 



|e, 0) - \l (e"«*f + l) |e, 1) + (e^- - l) \g, 1), 

(23) 

where U (0\g/e,0) = cos (§£) \g/e, 0) ±sin (| £) |e/ 5 ,0). 

The interference between evolution under W\ and W2 
is clearly seen in the amplitudes of |e,l) and \g, 1). 
For large squeezing, i.e. s = c, and qb — 0, we get 
|e) ® (JO) — -A|l)) which contains no entanglement. For 
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4> = 7r/2 we get |e, 0) + X\g, 1) which is entangled with 
entanglement strength on the order of the interaction 
strength A. This result is similar to that obtained for the 
tangle in III B2. For other atomic initial states similar 
calculations up to second order in A were performed 
leading to the results shown in III B2 and III B3. 

4b. Optical Block equation (OBE) with a noisy field 
An intuitive understanding of the error results from 
Sec. Ill B can be obtained from the Bloch sphere 
picture. To this end we derive the OBE for a one 
mode SL field, starting from the Hamiltonian (|SJ), H — 
ihg Ua + ct)a+ — {a) + a*)<7_] . As usual for SL, we de- 
fine the quadrature operators Aj. 2 — SX12 + £1,2 with 

a = i (sxt + isx 2 } 

a = ~ (x 1 +ix 2 ) ■ (24) 

Here 5Xi i2 are hermitian operators with zero average 
and Xi t 2 are real numbers representing the averages of 
the quadratures. We then define the vectors 

x = (x 2 ,xi,0) 

sx = (sx 2 ,5Xi,a) , (25) 

so that the Hamiltonian ([5]) becomes 

H = -h~(x + 5Xy%, (26) 

with a = (a xi a y ,a z ). The resulting Heisenberg equation 
for 3{t) is 

% = -lSxff-UxxS, (27) 
dt 2 2 y J 

with all the operators written in the Heisenberg picture. 
Rescaling the time to i = t/r Eq. (|2T|) becomes 

^ = -fi x S - ^(SX x S, (28) 

where f2 = (grx/2) is the dimensionless Rabi vector and 
C = gr. Recall that for a 7r pulse |a| = (ir/2gT) and 
therefore |f2| = 0(1). Our calculation also assumes A = 

cosh(r)C <C 1, i.e. even for matrix elements of SX which 
are of order cosh(r), (SX is still very small. We can 
then expand the operators a(t) = ao(i) + C<?i(i) + •■■ and 
5X(t) = 5X (t) + (5Xi(i) + to obtain to first order in 

c, 

^^-Qxff-l-CSXoxffo, (29) 



The last term includes the cross product of the Heisen- 
berg picture operators in the absence of atom-field inter- 
action, so that ao(t) is simply the solution of the OBE 
with the classical field ti, and describes evolution in the 
atomic part of the Hilbert space, and 5Xq — dX(i = 0) 
is in the photonic part of the Hilbert space. 

The average over atomic degree of freedom can be de- 
fined on any operator O to be (0)a = ^a[p"aO], where 
A denotes the atom system and pa is the reduced atomic 
density operator. Averaging equation (|29|) this way we 
get an OBE with a noisy field, 

a -fi x 0) A - ^8% x (ff ) A . (30) 

Note that (a) Ait) and SX are operators in the photonic 
space and are therefore fluctuating, whereas (t? )A(f) is 
just the c-number Bloch vector solution of an OBE driven 
by the average field Cl and it in fact equals to the average 
of ait) over both atomic and photonic degrees of freedom. 

The term with 5X$ is the one driven by the fluctuations of 
the field, which in our case are quadrature dependent SL 
fluctuations. For = 0, the uncertainties of the quadra- 
tures are A<SAi = e~ r < 1 and A<5A 2 = e r > 1, i.e. SXi 
has reduced fluctuations and 5X 2 has increased fluctua- 
tions relative to CL. For tj> = tt/2 we get the opposite. 

We can now present the field in the Bloch sphere as 
seen in Fig. 4. Taking real a (9 = 0) and squeezing 
phase (f> = 0, the Rabi vector is pointing to the y direction 
ix 2 = 0) and is amplitude squeezed (ASXi < 1). Now 
consider an initial ground state, i.e. (o) A is pointing to 
the —z direction. For a 7r pulse operation, the xy phase of 
the field is not relevant, and therefore phase noise of the 
field doesn't affect the operation. The only relevant noise 
is the amplitude noise which for <p = is reduced relative 
to CL case. We thus expect to get a lower error in this 
operation as compared with the CL operation. For <fi = 
it/2, the field has an increased amplitude noise (A£Xi > 
1), which results an increased error compared to the CL 
case. Now consider the y axis pointing initial state again 
for = 0. Since the Rabi vector is also pointing to the y 
direction the average "moment" obtained from the cross 
product is zero. Here only the 8X 2 "phase" noise, in 
the x direction, affects the qubit state. When this noise 
is increased (</> = 0) the error is larger than in the CL 
case, and when it is reduced {<p = 7r /2) the error gets 
lower. For an initial state pointing to the x direction, 
one can see that in general, both SX 2 and 8X\ affect 
the operation, and therefore squeezing in any quadrature 
typically increases the operation error. 

C. Entanglement and error probability calculations 
in the presence of the reservoir 

We now wish to examine how accounting for the 
spontaneous emission into the reservoir of empty modes 
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FIG. 4: The driving field in the optical Bloch equation repre- 
sentation. The solid arrow is the Rabi vector Q. proportional 
to x from Eq. (|25[l for X2=0 (6 = 0). The error bars are 
the variances of the field fluctuations ASX2 and ASXi in 
the x and y directions respectively. Here (j> = is taken, i.e. 
ASXi < A8X 2 . 



will modify the above results. To this end the full 
master equation, Eq. (fl3|) . has to be solved. Wc 
again assume A -C 1 and work in the truncated basis 
{|g,0), \g, 1), |e,0), |e, 1)}. We solve Eq. (fT3]) numerically 
in this basis using k — T/1000 and calculate the density 
matrix p(r) of the bipartite final state. We then use p{r) 
to evaluate the tangle and the error probability. 



1. Tangle calculation 

Within the truncated basis, the pulse mode is also a 
two-level system. We follow [2J, [25| where a procedure 
for the tangle calculation of a two-qubit mixed state is 
given. Figure [5] presents the tangle as a function of the 
pulse duration r. We assume squeezing of r = 1.16, 
which corresponds to 101og 10 (e 2r ) = lOdB pulse mode 
quadrature noise reduction. The interaction with the 
reservoir changes atom-pulse entanglement when Tt > 1 
and should generally decrease it for Tr >■ 1. Whenever 
the tangle is not suppressed by the SL operation, pulses 
shorter than O.ir -1 are well accounted by neglecting the 
reservoir, i.e. extending the initial linear part of the plots. 
When SL operation suppresses the entanglement, the in- 
teraction with the reservoir becomes more dominant, and 
the effect of the reservoir can only be neglected for pulses 
shorter than O.Oir^ 1 . 



2. Error probability calculation 

The error probability is calculated using Eq. (j2"Tj) . Fig- 
ure [S] shows the results for r = 0.1T -1 pulse duration. 



o 




- |g)+|e>, <|>=0 

- |g>+|e>, 1(1=11/2 
i|g>+|e>, *=0 
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FIG. 5: The tangle generated between the atom and the 
squeezed light pulse mode as a function of the pulse dura- 
tion scaled to atom free space decay time F _1 . The squeezing 
of the pulse mode is taken to be r = 1.16, which corresponds 
to 101og 10 (e 2r ) = lOdB quadrature noise reduction. As in 
Fig. 1, k = T/1000 is taken. The reservoir of the free space 
EM modes is included. Results for two atomic initial states 
are plotted. The effect of the reservoir is seen in the deviation 
from the initial linear part of the plots. This deviation be- 
gins as expected at times slightly shorter than r ~ T -1 , when 
entanglement is enhanced by SL (all cases here, apart from 
i\g) + |e), <j> — it/2). When entanglement is suppressed by SL, 
the effect of the reservoir is seen already at much smaller r. 



Let us now compare between Fig. [6] and Fig. 2, where 
the reservoir has been neglected. We focus on the case 
of the initial state 4= (i\g) + |e)). For r = 0, i.e. a CL 
pulse, the error of about 0.0244 seen in the left hand part 
of Fig. [6] is actually due to spontaneous emission to all 
the EM modes, as can be understood from the Mollow 
picture (see Ref. @). The corresponding error in Fig. 
2 is only 10~ 5 . This suggests that owing to its small 
solid angle (k = r/1000), the error due to coupling to 
the laser mode is overwhelmed by that due to coupling 
to the reservoir modes. When squeezing is introduced 
(r > 0), for <j> — 7r/2, the atom-pulse entanglement is 
suppressed as seen in Fig. 2. This suppression is hardly 
noticeable when the reservoir modes are included as one 
can see from Fig. [6] Numerically, for r — 3, we find that 
the error decrease is 0.95 x 10~ 5 . This matches what 
could be expected from Fig. 2 and demonstrates that 
the reservoir is the major cause of error in this case. On 
the other hand, for cp = 0, an increase of the error seen 
in Fig. 2 persists also in the presence of the reservoir, 
as can be seen in Fig. [SJ For very large squeezing this 
increase may become of the same order of magnitude as 
the error due to the spontaneous emission to all the EM 
modes. E.g. for r — 3 the increase seen in both Figs. 
2 and E] is ~ 4 x 10" 3 , about 16.3% of the 0.0244 error 
induced by spontaneous emission. 
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FIG. 6: The error probability of an operation with a squeezed 
light 7T pulse on an atom vs squeezing strength. The pulse 
duration is r = O.IT -1 and atomic decay into all the reservoir 
modes is considered. As in Fig. 1, k = r/1000 is taken. The 
left plot is the calculation for i\g) + |e) atom initial state and 
the right one is for \g) + |e). The dashed lines are for (j> — 
and the solid lines for <f> = 7r/2. 



IV. DISCUSSION 

The results of Sec. Ill show that the entanglement be- 
tween the pulse mode and the atom is indeed modified 
in a squeezing-phase-dependent way. However, the use of 
squeezed light would not reduce the average error com- 
pared with the CL case by a practicably useful amount. 
Despite its inapplicability for Quantum computing, the 
entanglement between SL and a single atom is of scientific 
interest, and would be interesting if observed experimen- 
tally. 

Contrary to CL, where atom-light entanglement is 
probably too small to be detected, SL with a large squeez- 
ing parameter r and a right choice of phases, will even- 
tually enhance SL-atom entanglement to experimentally 
detectable values. A direct measurement of entangle- 
ment would be through complete state tomography of 
the combined SL-atom state. The tangle would then be 
calculated from the reconstructed density matrix. The 
requirements for this experiment would be a squeezed 
light source with a sufficiently large squeezing parameter 
|28j , atomic state detection and photodetectors with suf- 
ficient efficiency. Alternatively, ignoring the light mode, 
the error in the final atomic state can be measured for 
various parameters, such as different squeezing phases, 
thus relaxing the need for an efficient photodetector. 

Examples for squeezed light sources include four 
wave mixing in a nonlinear medium, optical paramet- 
ric amplifiers, second harmonic generation and amplitude 
squeezed light from diode lasers among others [29(. All 
these sources can be pulsed, e.g. via electro-optic modu- 
lators. The initial light state we assumed for the control 
field was a single-mode squeezed state. Parametric down 
conversion (PDC) of short pulses could be considered as 
one example of a realization for such a light source. In 
[lil HH the multimode pulsed SL produced by PDC is an- 
alyzed and described by a discrete set of squeezed modes. 
In fact, for some conditions of phase matching and pump 
pulse duration, only one of the modes of the discrete set 



is squeezed, resulting in a pulse in a single-mode squeezed 
state. 

Even the most advanced single-mode squeezed light 
sources to-date cannot reach quadrature noise suppres- 
sion greater than lOdB (r = 1.16) [29]. This means 
that the maximal error increase and tangle expected and 
therefore the minimal required detection error, will be 
of the order of a few 10 -4 (see Fig. 1 and 2). Current 
state-of-the-art atomic state detection has an error which 
is of similar magnitude [32l . [33j . However the efficiency of 
even the most advanced photodetectors is far from sat- 
isfying such low error rates. It therefore seems that a 
direct tangle measurement is currently unfeasible, how- 
ever, an indirect indication for entanglement, detected 
through gate error enhancement, is within experimental 
reach. 

One important limitation for observing the error mod- 
ification effects in free space is the spacial overlap be- 
tween the atom and the paraxial pulse mode. Similar 
limitations make up a major difficulty in observing the 
effect of the inhibition of the atom coherences decay in 
a squeezed vacuum reservoir, found by Gardiner [l5j . In 
this context, other schemes, involving atom-squeezed vac- 
uum interaction in a cavity, were considered [34j | . Recent 
theoretical analysis and experiment however have shown 
a possibility of strong coupling of a laser mode to a single 
atom in free space [35l l36l |37I |. The use of such beams, 
where the overlap between the incoming beam and the 
atomic dipole radiation pattern is large (re ~ T), may 
greatly enhance the effects discussed above. 

In many experiments, a pair of ground, e.g. hyperfinc, 
states of an atomic system is used as a qubit. Although 
the energy separation between ground states is typically 
in the radio-frequency range, quantum gates are often 
driven with a laser beam pair, off-resonance with a third 
optically excited level, in a Raman configuration. The 
limitations to hyperfine qubit coherence and gate fidelity 
due to the quantum nature of the Raman laser fields 
was investigated both experimentally and theoretically 
[H, [39[ . In this case too, it is interesting to ask 
whether the use of, either independently or relatively, 
squeezed light fields would influence the gate error. An 
experiment, similar to that described above, but with 
a ground state qubit could benefit from the ground 
state qubit intrinsic longevity. In fact, a gate averaged 
error of few 10~ 3 , dominated by laser classical noise and 
photon scattering, was recently measured on a single 
trapped-ion hyperfine qubit and a Raman laser gate [40( . 

To conclude, we have calculated the entanglement be- 
tween a short pulse of resonant squeezed light and a two- 
level atom in free space during a qubit ia y gate (a tt 
pulse) and the resulting gate error. Several differences 
arise when comparing to coherent control with CL. As 
could be expected for squeezed light, three phases deter- 
mine the error magnitude - the squeezing phase, the light 
phase and the phase of the initial qubit superposition. In 
comparison with the coherent light case, the error and 
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the entanglement can be either enhanced or suppressed 
depending on the above phase relations. These entan- 
glement enhancement and suppression effects naturally 
become stronger with squeezing and can get to a few 
orders of magnitude. For quantum information process- 
ing purposes the relevant error is that averaged over all 
possible initial states. We have shown that although in 
principal the average error can be reduced by using SL 
within a specific range of parameters, the resulting min- 
imum error however would be of the same order of mag- 
nitude as that with CL. In fact, even when not averaged, 
but when taking into account the free space environment, 
the error is still mostly dominated by the atomic decay 
to the EM reservoir owing to a very small atom-pulse 
spacial overlap, so that any reduction is not likely do be 
detected experimentally. On the other hand, the effect 
of enhanced atom-pulse entanglement is of scientific in- 
terest and might even prove useful as a resource in some 
quantum communication schemes. A measurement of the 
increased error would be an indirect indication for the en- 
hancement of entanglement by squeezing. 
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the modes centered around some carrier wave vector in 
the z direction k a = k^z: 



{A*,A£,A*} 

a; = (fco-fcf,fco- 



p 

A" 

p 



(A.2) 



where the bandwidths for (x,y,z) directions are narrow, 
i.e. 



k^ ko 
i = x,y,z. 



(A.3) 



The nonparaxial modes are simply the rest of the modes 
k\x not included in A p and we denote them with indices 
{q}. The paraxial subspace A p can be approximately 
spanned by the modes which are the solutions of the 
paraxial wave equations in free space with carrier fre- 
quencies uik = {ko + k)c, TEM k °+ k . These modes are 
denoted by the functions U nm k{r) — u nm k(r)e l ^ k ° +k ' z . 
Their Fourier coefficients are U n mk(k) so that their mode 
operators are 



We appreciate useful discussions with Ilya Averbukh, 
Eran Ginossar and Barak Dayan. RO acknowledges sup- 
port from the ISF Morasha program, the Chais family 
Fellow program and the Minerva foundation. 



APPENDIX: QUANTIZED PARAXIAL 
APPROXIMATION 



Here we would like to derive the model Hamiltonian 
of the atom-pulse mode interaction in free space seen in 
Eq. |J5J| . We start from the usual atom- field Hamiltonian 
in the dipole approximation 



H 

Hf 

H A 
Haf 



H F + H A + H AF 

1 



hw a a z 



E 

k fx 



%\ / at. ar 



ik-r a 



h.C. (<7_|_ + <7_) . 



(A.l) 



where kfj, are the indices of a Fourier mode and its polar- 
ization respectively, and r a is the position of the atom. 
The dipole matrix element d is assumed to be real and 
df. = eV • d, where eg is the polarization vector with 
polarization /i and is perpendicular to the propagation di- 
rection k. We first divide the sums on Fourier modes into 
the paraxial subspace modes, A p , and the nonparaxial 
subspace modes. The paraxial subspace modes include 



keA P 



k fi 



nrnk 



(k)a. 



nmkjfj, ■ 



(A.4) 



Note that the Fourier modes here, ag , are only those 
included in the paraxial subspace. Using (|A.4|) it is easy 
to show within the paraxial approximation that 



ik-r* a 



v 



wo E d ^ 

nmk,fi 

(A.5) 

where u>o = k c. We denote all the {nm ^ 00, k, n} and 
{nm — 00,k,u — 2} modes by {«}, where {fi} = {1,2} 
are the two polarizations. The rest of the modes ,{nm = 
00, = 1}, have Gaussian transverse profile and are 
denoted {k}. We thus get 



H A i 



«>E 



9 

E 



9q a q 



9* q a\) (<7+ 



iH 



{9iai-gta?J (jx + +<7_) + 



k 



^0 

2e h 



dU ok(r"a)ak - h.c (6+ + <r_) 



(A.6) 



where d = d^ = \ and g q are the atom-mode dipole 
couplings for the modes {i,q}- The dipole couplings for 
the {k} modes are seen explicitly in (|A.6j) . 



9k 



2e Q h 



dU 00k (f a ) 



^0 

2e ?f 



du 



mk (ray (ko+k)Za - 

(A.7) 
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We note that Moofc(fa) is just the overlap of the fc mode 
with the atom. This overlap includes some normalization 
of a quantization length in the beam direction, L, which 
should eventually be taken to be infinity, and which is 
multiplied by some effective area of the mode, A, eval- 
uated at the atom's position, i.e. tioofc(^a) = 1/yAL. 
Considering the narrow band of fc values around ko, A 
is approximately independent of k. The couplings g k are 
then written 



9k 



1 



2e hA VI 



a i(k +k)z a 



(A. 



As for Hp, since all the modes discussed here, {q, i,k}, 
are eigenmodes of free space in the paraxial approxima- 
tion, they have well-defined single-particle energies f)io q , 
hiOi and fkuk respectively. The total Hamiltonian from 
Eq. (|A.1[) thus becomes 



H 

Hf 

H A 
Haf 



H F + H A + H AF 

^ fiw q a q a q + ^2 fruia\ai + ^ frw k a\a k 

q i fc 

ih {so & j - 9*jafj (<r+ + o--) • (A.9) 

j={q,i,k} 



We would now like to describe a pulse mode with a 
Gaussian transverse profile composed of Gaussian beams 
of different frequencies around the carrier frequency u>q. 
The pulse mode is then spanned by the {k} modes. We 
can choose a new orthonormal basis {/} which spans the 
same subspace as that of the {fc} basis and that includes 
a pulse mode fo as one of its members. Both {fc} and {/} 
modes share a similar Gaussian xy profile, we therefore 
focus only on their z dependence. We write the relation 
between the k and / mode functions as 



4>, 



k 



„ikz 



s/L' 



(A.10) 



where e^ ko+k ^ z / y/Z is a fc mode while $/(z) is an / mode. 
An / mode is thus defined by its envelope ipf(z). In a 



similar way, the transformation for the field operators is 



(A.ll) 



/ 



The pulse mode is then characterized by the operator 
a/ , whereas the Hilbert space of the two-level atom is 
spanned by the Pauli spin matrices. The atom and the 
pulse mode make for the relevant bipartite system. The 
rest of the EM modes, S/^/ , Sj and a qi interact with the 
bipartite system via the atom and make up a reservoir. 

We now focus on the part of the Hamiltonian that in- 
cludes only the atom and the Gaussian paraxial modes 
{fc}, and their interaction 

H p = hhj k a\a k + ^friv a a z + 
fe 

ihJ2{g k a k ~h.c.) (<7++<7_). (A.12) 



Moving to the interaction picture in the rotating wave 
approximation, assuming no detuning (ui a = luq) and us- 
ing equations (|A.8[ lA.lOi lA.llj) . we get 



Hf = ih^2[ g f (i)a f a + - h.c] , 
/ 



(A.13) 



with 



9 fit) 



2e HA 



e lkoZ ^ f (z a -ct). 



(A.14) 



The system Hamiltonian of the relevant bipartite system 
comprised of the atom and the pulse mode is then 



H S = Hi 



fo-\-atom 



ihg fo (t) [ a fo &+-h.c], (A.15) 



where we take the envelope of the pulse, <ff (z), to be real 
and set z a = 0. Denoting all the rest the EM modes apart 
from fo with indices {r} we get the total Hamiltonian in 
the interaction picture in Eq. @. 
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